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ABSTRACT 

We have studied the properties of the stellar populations in the field of the NGC 346 cluster in the Small 
Magellanic Cloud, using a novel self-consistent method that allows us to reliably identify pre-main sequence 
(PMS) objects actively undergoing mass accretion, regardless of their age. The method does not require spec- 
troscopy and combines broad-band V and / photometry with narrow -band Ha imaging to identify all stars with 
excess Ha emission and derive the accretion luminosity La^c and mass accretion rate M^^c for all of them. The 
application of this method to existing HST/ACS photometry of the NGC 346 field has allowed us to identify 
and study 680 bona-fide PMS stars with masses from ~ 0.4 Mq to ^ 4 M0 and ages in the range from ~ 1 Myr 
to ~ 30 Myr Previous investigations of this region, based on the same data, had identified young (~ 3 Myr old) 
candidate PMS stars on the basis of their broad-band colours. In this study we show that there are at least two, 
almost equally numerous, young populations with distinct ages of respectively ^ 1 and ^ 20 Myr We provide 
for all of them accurate physical parameters. 

We take advantage of the unprecedented size of our PMS sample and of its spread in mass and age to 
study the evolution of the mass accretion rate as a function of stellar parameters. We find that, regardless 
of stellar mass, the mass accretion rate decreases with roughly the square root of the age, or about three 
times slower than predicted by current models of viscous disc evolution, and that more massive stars have 
systematically higher mass accretion rate in proportion to their mass. A multivariate linear regression fit reveals 
that logMacc — -0.61ogf + log m + c, where t is the age of the star, m its mass and c a quantity that is higher at 
lower metallicity. This result is consistent with measurements of the mass accretion rate in the 30 Dor region 
and in the Milky Way and suggests that longer duration for mass accretion could be related to lower metallicity. 
The high mass accretion rates that we find suggest that a considerable amount of mass is accreted during the 
PMS phase, of order ^ 0.2 Mq or possibly ^ 20 % of the final mass for stars with mass m < 1 Mq if their discs 
are eroded by 20 Myr, i.e. before they reach the main sequence. Therefore, PMS evolutionary models that do 
not account for this effect will systematically underestimate the true age when compared with the observations. 

Subject headings: accretion, accretion disks - stars: formation - stars: pre- main-sequence - Magellanic Clouds 



1. INTRODUCTION 

NGC 346 is presently the region of most intense star for- 
mation in the entire Small Magellanic Cloud (SMC). It con- 
tains over 30 O type stars (Massey, Parker & Garmany 1989; 
Evans et al. 2006) that ionise the N66 nebula, the largest 
HII region in the SMC (Henize 1956). The presence of these 
very young (^1-3 Myr) massive stars in NGC 346 has long 
been established with ground-based spectroscopy (e.g. Wal- 
born 1978; Walbom & Blades 1986; Niemela et al. 1986; 
Massey et al. 1989), but more recently also very young stars 
of low mass have been detected. Observations with the Hub- 
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ble Space Telescope and Spitzer Space Telescope have re- 
vealed respectively a multitude of pre-main sequence (PMS) 
star candidates down to the subsolar mass (Nota et al. 2006; 
Sabbi et al. 2007; Gouliermis et al. 2007; Hennekemper et al. 
2008) and a large number of candidate young stellar objects 
(YSO) with masses as low as 1 .5 M© (Bolatto et al. 2007; 
Simon et al. 2007). 

What makes NGC 346 and the surrounding regions particu- 
larly interesting is that they allow us to study the properties of 
star formation in a nearby galaxy with a metallicity similar to 
those in place in the high redshift universe at z > 2. The dis- 
tance modulus for the SMC is 18.92 ± 0.03 corresponding to 

- 61 kpc (Hilditch et al. 2005; Keller & Wood 2006), whereas 
the currently accepted values for its metallicity range from 

- 1/5 to - 1/8Zq (see Russell & Dopita 1992; Rolleston 
et al. 1999; Lee et al. 2005; Perez-Montero & Diaz 2005). 
Moreover, besides this recent burst, previous star formation 
episodes have been detected in this area. Heap et al. (2006) 
and Mokiem et al. (2006) reported the presence of massive 
stars with an age of ^ 5 Myr, while Massey et al. (1989) 
discovered five red supergiants and two B type supergiants 
that form a spatially distinct subgroup located about 2. '6 SW 
of the centre of NGC 346 and that have an estimated age of 
^ 12 Myr More recently, Sabbi et al. (2007) detected a small 
star cluster, located about 2.'1 or '--^ 37 pc NE of the centre 
whose colour-magnitude diagram (CMD) is compatible with 
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an age of 15±2.5Myr. Finally, older populations are also 
present: Sabbi et al. (2007) find evidence for star formation 
in this region dating as far back as ^ 10 Gyr with a moderate 
enhancement ^ 150 Myr ago. In addition, an intermediate- 
age cluster, BS 90, is also present in the field, with an age of 
4.3 ±0.1 Gyr 

It is not easy to establish whether there is any relationship 
between the very old generations in this field and the more 
recent star formation episodes. However, it is interesting to 
investigate whether the ^ 12 Myr old massive stars of Massey 
etal. (1989) and the 15 ±2.5 Myr old small cluster of Sabbi et 
al. (2007), although spatially unrelated to one another, could 
be the precursors of the 30 young massive O type stars and 
of the recently detected PMS stars and YSO in this region. It 
is conceivable that a large number of these ~ 15 Myr old ob- 
jects are actually present in NGC 346, many more than those 
known so far, but that they could belong to a more diffuse 
population spread throughout the region because of apprecia- 
ble velocity dispersion motions and, therefore, might not be 
easy to identify. If the members of such a 15 Myr old pop- 
ulation could be detected, including those of low mass, their 
spatial distribution and physical properties would allow us to 
understand whether and how they may have triggered the most 
recent bursts, as theories of sequential star formation suggest 
(e.g. Elmegreen & Lada 1977). 

The task of detecting a ^ 15 Myr old population in a region 
heavily contaminated by older field stars and younger objects 
might appear daunting. Hennekemper et al. (2008) attempted 
to study the presence of an age spread in the PMS population 
of NGC 346 through the analysis of the CMD. They showed 
that the observations are in principle compatible with an upper 
limit of ^ 10 Myr to the age, but they concluded that differen- 
tial reddening and unresolved binaries could be at the origin 
of the observed colour spread. 

Actually, if these objects have not yet completed their pre- 
main sequence (PMS) phase, it should be possible to detect 
them through the distinctive excess emission features in their 
spectra that originate from the accretion process, particularly 
in the UV continuum and in recombination lines such as Ha, 
Pa/3 and Br7 (e.g. Calvet et al. 2000). Although in nearby star 
forming regions the fraction of PMS stars with discs appears 
to decline quite rapidly during the first 10 Myr (e.g. Haisch, 
Lada & Lada 2001; Fedele et al. 2010), there is evidence 
that sustained mass accretion is present also in objects older 
than this age, particularly above 1 Mq . This is the case for 
instance of Tr 37 where Sicilia-Aguilar et al. (2006) find G- 
type dwarfs still accreting at ages in excess of 10 Myr In 
the more distant (^ 7kpc) and massive Galactic star form- 
ing region NGC 3603 Beccari et al. (2010) recently found 
a conspicuous number of PMS stars older than 10 Myr still 
undergoing mass accretion. In even more distant and dense 
stellar fields, where spectroscopy of individual stars is lim- 
ited to the brightest members, Romaniello (1998), Panagia 
et al. (2000) and Romaniello, Robbeito & Panagia (2004) 
have shown that objects undergoing active mass accretion can 
be efficiently detected with accurate multi-colour photometry 
and have identified in the regions around SN 1987 A about 500 
PMS stars with Ha equivalent width in excess of 8 A and with 
an age of 12 ±2 Myr More recently, De Marchi, Panagia & 
Romaniello (2010, hereafter Paper I) showed that through a 
suitable combination of broad- and narrow-band photometry 
it is also possible to derive the mass accretion rate of these 
objects, with an accuracy comparable to that allowed by spec- 
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Figure 1. The obsei'ved colour-magnitude diagram of a field of ~ 200" on 
a side containing the cluster NGC 346 in the Small Magellanic Cloud, as 
observed with the Advanced Camera for Surveys on boai'd the Hubble Space 
Telescope by Nota et al. (2006). The diagram, based on the photometric study 
of Sabbi et al. (2007), contains and a total of about 53 550 stars detected 
at the > 3 (T level in V and /, of which 1 8 764 have a combined magnitude 
unceitainty in the V, I and Ha bands of less than 0. 1 mag (see Equation[T). 

troscopy. 

In this work, we build on the method developed in Paper I 
to securely identify bona-fide PMS objects in dense stellar 
fields and apply it to the high-quality HST photometry of the 
NGC 346 region (Sabbi et al. 2007). Our goal is to mea- 
sure the physical properties of these objects and study how 
they depend on age, in order to understand how star forma- 
tion has proceeded in this area over the past ~ 30 Myr The 
paper is organised as follows: in Section|2]we briefly describe 
the observations and address how we correct for differential 
reddening, whereas Section|3]is devoted to the identification 
of PMS stars via their Ha excess emission and on the deter- 
mination of their Ha luminosity. In SectionH] we derive the 
other important physical parameters for these objects, includ- 
ing age, mass and mass accretion rate, while in Section|5]we 
study how the latter evolves in time, how it could affect stel- 
lar evolution, and we compare the results to those for similar 
objects in the Milky Way and Large Magellanic Cloud. A 
summary of the most important conclusions of the paper is 
offered in Section|6] 

2. THE OBSERVATIONS 

The observations used in this work were collected with the 
Wide Field Channel of the Advanced Camera for Surveys 
(Ford et al. 2003) on board the Hubble Space Telescope in 
2004 July as part of proposal GO-10248 (with AntoneUa Nota 
as principal investigator). The data and their reduction were 
extensively described in Sabbi et al. (2007), to which we refer 
the reader for details. In summary, the observations of inter- 
est here cover an area of 200" on a side centered on the 
cluster NGC 346 in the Small Magellanic Cloud (SMC), at 
RA= 00''59'"5.2^ DEC=-72°10'28" (J2000), itself immersed 
in N 66, the largest and brightest HII region in the SMC. 

The data comprise both deep and short exposures in the 
broad bands F555W (hereafter V) and F814W (hereafter /), 
as well as observations in the 80 A wide band F658N (here- 
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after Ha), which were not described in Sabbi et al. (2007). 
The total combined exposure times in each band are, respec- 
tively, 4 148 s, 4 124 s and 1 542 s. The data reduction in the 
y and / bands was carried out by Sabbi et al. (2007) using a 
combination of PSF fitting and aperture photometry and their 
final catalogue, in the central 200" x 200" of interest here, in- 
cludes about 53 550 stars reliably detected (at the > 3 cr level) 
in y and /. As regards the Ha band, we have adopted the 
same combination of PSF fitting and aperture photometry to 
search for and measure all objects with a signal in the central 
pixel exceeding the local background level by at least 3 times 
the standard deviation of the latter (i.e. > 3cr detection). Of 
all the stars detected in the V and / bands, about 19 800 satisfy 
this condition also in the Ha band. 

We show in Figure[T] the resulting colour-magnitude dia- 
gram (CMD) of all the stars in the selected field as light grey 
dots and with darker symbols stars whose mean error Si in the 
three bands does not exceed 0.1 mag, where 

and (^555, (5658 and Jgn are the photometric uncertainties in 
each individual band. A total of 18764 stars satisfy this con- 
dition, with the most stringent constraint being set by the un- 
certainty in the Ha band. 

The CMD, as already discussed by Nota et al. (2006), 
GouUermis et al. (2006; 2007) and Sabbi et al. (2007), reveals 
the presence in this field of both young and old populations. 
The former includes the massive stars on the pronounced up- 
per main sequence (MS) and for the candidate PMS objects 
to the right of and brighter than the lower MS, while the latter 
accounts for the lower MS and for the sub-giant and red giant 
branches (SGB, RGB). In general, however, it is not possible 
to classify an object as being young or old just on the basis of 
its position on the CMD, not only because there are regions 
such as the lower MS that are common to both populations, 
but also because there is evidence for variable reddening in 
the field, as Sabbi et al. (2007) already pointed out and as one 
can see from the conspicuous broadening of the upper MS in 
Figure[T] 

Differential reddening makes it particularly difficult to un- 
derstand the true nature of the large group of objects below the 
SGB and redder than the MS. While their colours and mag- 
nitudes, as originally pointed out by Nota et al. (2006), are 
compatible with those of young PMS stars, depending on the 
amount of reddening and on the specific extinction law for this 
field, some of them may in reality be reddened MS or SGB 
stars. Furthermore, even if all these objects were PMS stars, 
without accounting for differential reddening one could not 
properly assign to them a mass and age through the compar- 
ison with PMS evolutionary tracks (see e.g. Hennekemper et 
al. 2008). It is therefore important to address the issue of dif- 
ferential reddening and to correct the photometry accordingly, 
to the extent possible. We do so in the following section. 

2.1. Reddening correction 

Sabbi et al. (2007) concluded that the young population in 
this field, corresponding to the upper MS, is compatible with 
an age of 4 ± 1 Myr Therefore, one can determine the red- 
dening towards each of those stars if one assumes that they 
are all MS objects and one looks at their colour and magni- 
tude displacement with respect to a theoretical isochrone for 



an age of 4 Myr, a metallicity of Z = 0.002 and a distance 
modulus of 18.92 ±0.03 (see Introduction). Hennekemper et 
al. (2008) followed a similar approach to derive an extinc- 
tion map for this region, albeit with a set of isochrones for a 
higher metallicity (Z = 0.004). We note here that although the 
metallicity that we adopt, Z = 0.002, is at the lower end of the 
currently accepted values for the SMC (see Introduction), it 
appears better suited to describe the properties of young PMS 
stars, as we will show in Section|4] 

As for the extinction law in this field, we decided to use the 
one derived by Scuderi et al. (1996) for the field of SN 1987A 
in the Large Magellanic Cloud, rather than the Galactic law 
as adopted by Hennekemper et al. (2008), in light of the more 
similar metallicity. The total extinction values for the spe- 
cific HST filters used here were derived by averaging over the 
entire bandpass the dust attenuation (10"" '*'*^), as explained 
in Romaniello et al. (2002). The solid lines cutting through 
the upper MS in Figure|2^ correspond to the reddening vectors 
drawn from a 4 Myr old isochrone for Z = 0.002 (Marigo et al. 
2008), in steps of 0.5 mag for stars brighter than yobs — 18.5. 
The length of the lines corresponds to Ay = 2. A total of 
about 600 stars selected in this way and bluer than V -I = 0.3 
are marked with thicker symbols in Figure|2^. The corre- 
sponding reddening distribution ranges from E(V -I) = 0.1 
to E(V-I) = 0.25 (respectively 17 % and 83 % limits corre- 
sponding to Ay values of 0.27 and 0.67), with a median value 
of E{V-I) = 0.16 or Ay = 0.43, indicating that a considerable 
spread in reddening is present in this field. 

We have used the reddening values towards each of the se- 
lected upper MS stars to derive a reddening correction for all 
other objects in their vicinities. The analysis of the Ha image 
indicates that the typical scale of the brightness fluctuations is 
of order 5" radius. We have therefore assumed that stars pro- 
jected within this region have similar extinction, even though 
we have no information on the depth of their distribution nor 
on the uniformity of the absorbing material. In order to vali- 
date the procedure, we first applied the reddening correction 
to the upper MS stars themselves, using for each object the 
average of the extinction values of at most 20 other upper 
MS stars within 5" of their location. We took the straight 
average without weighting them for the distance from the ob- 
ject, since we have no information on the distribution of the 
stars along the line of sight. We found that the average and 
standard deviation of the reddening distribution derived with 
the nearest neighbours are in excellent agreement with those 
measured directly, namely < Ay >= 0.51 ±0.19. This has 
convinced us that the average reddening value measured in a 
small area surrounding each star can be used to correct for 
extinction. Objects with no reference upper MS stars within 
5" were assigned the average reddening value of all reference 
stars (Ay = 0.51). 

The upper MS corrected in this way appears considerably 
tighter in the CMD than the raw data (see FigurelJt) and this 
convinced us that the method, albeit statistical in nature, can 
be applied to the entire data set. We also experimented with a 
different number of neighbours, namely 5, 10, 20 and 50. We 
found that, by using the 20 nearest upper MS stars, the red 
clump (RC) feature in the red giant branch of old ('^1 Gyr) 
in the CMD at y ~ 19 appeared the tightest. Since our goal 
is to correct for reddening all candidate PMS stars, which are 
typically red, this choice seemed the most appropriate. The 
de-reddened CMD obtained in this way is shown in FigurellJ). 
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Figure 2. Panel a): Observed CMD as in Figure[T] with overplotted tlie isoclirone of Marigo et al. (2008) for an age of 4Myr, metallicity Z = 0.002 and distance 
modulus (m—M)Q = 18.92 ± 0.03. The slanted thin solid lines show the reddening vectors, for Ay = 2, according to the extinction law of Scuderi et al. (1996). 
They are drawn every 0.5 mag for stars brighter than V^bs ~ 18.5 and are used to identify likely stars in the upper MS that have been displaced by reddening 
(shown as thick dots). Panel b): the same CMD after correction for differential reddening, using as a reference for each object the average extinction value of the 
twenty nearest upper MS stars. Note that both the upper MS and RC are tight. 

3. IDENTIFYING THE PRE-MAIN SEQUENCE STARS 

One of the characteristic signatures of the accretion process 
on to PMS stars is the presence of excess emission in Ha (see 
e.g. Calvet et al. 2000). Recently, we showed in Paper! that 
PMS stars undergoing active accretion can be reliably identi- 
fied, and their Ha luminosity measured, by using a combina- 
tion of broad-band (V,I) and narrow -band (Ha) photometry, 
with no need for spectroscopy. This way of identifying PMS 
stars is more accurate and reliable than the simple classifica- 
tion based on the position of the objects in the Hertzsprung- 
Russell diagram, i.e. stars placed well above the main se- 
quence (e.g. Gilmozzi et al. 1994; Hunter et al. 1995; see in 
particular Nota et al. 2006 and Gouliermis et al. 2007 for the 
specific case of NGC 346). The true advantages of the method 
presented in Paper I are that: (/) it allows us to discriminate 
between bona-fide PMS stars and interlopers that occupy the 
same region of the CMD; it provides a secure detection 
of relatively old PMS stars, already close to their MS, whose 
colours and magnitudes would otherwise be indistinguishable 
from those of MS stars; by measuring the Ha luminosity 
L(Ha) of these PMS stars we can obtain their mass accretion 
rates Mace from L(Ha) and the stellar parameters (mass and 
radius); (iv) finally, this approach is more practical, efficient 
and economical than slitless spectroscopy for studying PMS 
stars in dense star forming regions, although it cannot provide 
constraints for magnetospheric accretion models that require 
information on the profile of the Ha emission line (e.g. Muze- 
rolle et al. 2001). 

As explained in detail in Paper I, the method hinges on the 
simple consideration that, at any time, the largest majority of 
stars of a given effective temperature Teff in a stellar popula- 
tion will have no excess Ha emission. Therefore, for stars 
of that effective temperature, the median value of any colour 
index involving Ha (e.g. V -Ha) effectively defines a spec- 
tral reference template with respect to which the Ha colour 
excess should be sought. 
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Figure 3. Colour-colour diagram of the selected 1 8 764 stars with 53 < 
0.1 mag. The dashed hne represents the running median V -Ha colour, ob- 
tained with a box-car size of 100 points, whereas the thin solid line shows 
the model atmospheres of Bessell et al. (1998). The arrow in the figure cor- 
responds to the reddening vector for Ay = 2.7 or E{V — I) = 1, taken as an 
example, for the specific bands used in this study. A total of 79 1 objects with 
a V -Ha excess larger than 4(T are indicated with large dots (in red in the 
onhne version). 

In Figure[3]we show with small dots the V -Ha colour as 
a function of V - / for all 1 8 764 stars that have a combined 
photometric uncertainty ^3 < 0.1 mag (see Equation[T]). The 
value of ^3 is dominated by the uncertainty on the Ha mag- 
nitude, while the median value of the uncertainty in the other 
two bands is Sy = 0.006 and 6f = 0.008, respectively. 

The thick dashed line in Figure[3] represents the median 
V - Ha colour obtained as the running median with a box- 
car size of 100 points. The thin solid line shows the colours 
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in these filters for the model atmospheres of Bessell et al. 
(1998). As discussed in Paper I, the < 0.1 mag discrepancy 
between models and observations around V- / ~ 0.5 is most 
likely due to the coarse spectral sampling of the Ha line in 
the models. Where the density of observed points decreases 
considerably, at V - / > 1 .2, we reduced the box-car size to 
10 points and the result is shown by the thick dot-dashed line 
(which is extrapolated for V -I > 2). A total of 891 objects 
have a y - Ha index exceeding that of the reference template 
(i.e. the photospheric colour) at the same V-I colour by more 
than four times the uncertainty on their V- //a values. 

The arrow in the figure corresponds to the reddening vec- 
tor for Ay = 2.7 or E{V-I) = 1 for the specific bands used in 
this study, derived from the extinction law of Scuderi et al. 
(1996). The fact that the arrow is almost parallel to the refer- 
ence template implies that even relatively large uncertainties 
in the reddening correction would not significantly affect the 
selection of stars with Ha excess. On the other hand, it is pos- 
sible that for some objects surrounded by non homogeneous 
nebular emission the Ha excess that we measure is spurious. 
We explain in Section lTTI how we have removed these objects 
from our list. 

3.1. Nebular contamination 

The effects of nebular contamination on the detection of 
stars with Ha excess emission are addressed in detail in Pa- 
per I (and in Beccari et al. 2010), to which we refer the reader. 
We show there that gas emission acts like a "pedestal," or an 
additional source of background. As such, it does not affect 
our photometry because it increases the measured flux in the 
same way both inside the central aperture and inside the back- 
ground annulus used for the determination of the sky level. 
Therefore the gas emission is automatically removed in the 
sky subtraction procedure. If the diffuse emission is strong 
compared to the intensity of the objects, the uncertainty on 
the flux of the star will be larger, but this is fully taken into 
account in our photometry. 

A possible complication, not discussed in Paper I, arises 
when the nebular emission is not homogeneous across the 
field. For example, when a filamentary structure in the neb- 
ular gas overlaps with a star, the magnitude of the star can 
be affected in different ways. If the filament is wider than 
the diameter of the annulus that we use to measure the local 
background, this only results in a higher photometric uncer- 
tainty, because we fall in the case explained above for diffuse 
emission. The situation is more complex when the filament 
only partly covers the background annulus (which in our case 
extends from 3 to 5 pixel radius around the object) and/or the 
central star (for which an aperture of 2 pixel radius is used), 
since in those cases the signal from the star and/or the back- 
ground could be incorrectly determined. 

We, therefore, inspected the Ha image around the location 
of all stars with an apparent Ha excess. In order to enhance 
the contrast between neighbouring pixels in the image, we ap- 
plied an unsharp-masking filter (e.g. Malin 1977) to the Ha 
frame, with a smoothing radius of 10 pixel meant to subtract 
the low-frequency signal, i.e. structures in the image with a 
characteristic scale larger than the size of our background an- 
nulus. Unsharp masking makes structures and filaments in 
the nebular gas easier to see than in the direct image. In this 
way we compared the structures and filaments with the po- 
sitions of the stars with Halpha excess emission (see Beccari 
et al. 2010 for an application of this method to the case of 
NGC 3603). We adopted a very conservative approach and 



excluded from our PMS candidate list all stars whose photo- 
metric background annulus is partially covered by a gas fil- 
ament. In total, about 70 stars were excluded in this way, 
or less than 10 % of the original population (an additional 30 
stars were excluded because they are closer than 10 pixels to 
the image edges in the Halpha frame, although their magni- 
tudes are well defined in the other bands). 

A total of 791 objects are left after this conservative selec- 
tion and they are indicated with large dots in Figure|3](shown 
in red in the on-line version) and, in light of their V - Ha ex- 
cess at the 4(7 level, they are excellent PMS candidates. As 
we shall see later, most of them are indeed bona-fide PMS 
stars. 

3.2. Ha luminosity 

Following the detailed procedure outlined in Paper I, we can 
derive the Ha emission line luminosity L{Ha) for these stars 
from the AHa colour corresponding to the excess emission, 
namely: 

AHa = (V- Ha)°^' - (V - Haf''^ (2) 

where the superscript obs refers to the observations and ref 
to the reference template. L(Ha) is then obtained from AHa, 
from the photometric zero point and absolute sensitivity of the 
instrumental set-up and from the distance to the sources. The 
assumed distance modulus is (m-M)o = 18.92 (see Introduc- 
tion), whereas the photometric properties of the instrument 
are as listed in the ACS Instrument Handbook (Maybhate et 
al. 2010). The derived median Ha luminosity of the 791 
objects with Ha excess at the Aa level is 2.7 x lO-'' erg/s or 
^ 10"^ L0. This value is similar to that of a group of 189 stars 
with Ha excess emission in the field of SN 1987 A that were 
studied in Paper I with the same method applied to HST Wide 
Field Planetary Camera 2 data. The uncertainty on L(Ha) is 
typically ^ 12 % and is completely dominated by the statisti- 
cal uncertainty on the Ha photometry (10%), while the dis- 
tance accounts for 5 % (Hilditch et al. 2005; Keller & Wood 
2006) and the absolute sensitivity of the instrumental setup 
for another 3 %. 

Owing to the rather wide passband of the F658N filter, 
AHa includes a small contribution from the [Nil] emission 
features at 6548 A and 6584 A. Even though such a contribu- 
tion is small (see Paper I, where the effect is explained in de- 
tail), it is a good practice to correct for it as it is a systematic 
effect. The adopted average correction factor to the intensity 
for the ACS F658N filter, based on the spectra of classical T 
Tauri stars is 0.979. 

3.3. Ha equivalent width 

Besides the luminosity of the Ha emission line, also its 
equivalent width Weq(//a) can be derived from our photome- 
try. To this aim, as shown in Paper I, we used the model atmo- 
spheres of Bessell et al. (1998) to derive the magnitude Ha'' 
corresponding to the sole continuum in the Ha band from the 
V and / magnitudes. The corresponding relationships, prop- 
erly validated by comparison with spectro-photometric mea- 
surements, are listed in the Appendix of Paper I and allow us 
to derive Weq(i/a) from the Ha magnitude via the equation: 

W^c^iHa) = RW X [1 - io-"-4x(»"-^^"')] (3) 

where RW is the rectangular width of the filter (similar in def- 
inition to the equivalent width of a line), which depends on 
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Figure 4. Ha equivalent width of tlie stars in the field of NGC 346, as a 
function of their V -I colour. Of the 791 stars with Ha excess at the 4 cr 
level (circles and squares), a total of 694 have Wcq(Ha) < —20 A (or < -50 A 
for stars hotter than 10 000 K). As such, we consider them as bona-fide PMS 
stars and indicate them with thick dots (shown in red in the online version). 
As customary, we use negative equivalent widths for emission lines. 

the characteristics of the fiher and corresponds to 74.96 A for 
the F658N fiher used here. The equivalent widths obtained 
in this way are shown in Figured as a function of the V -I 
colour The 791 stars with a 4a Ha excess are indicated with 
thicker symbols (circles and squares). We take as bona-fide 
PMS stars all those with Weq(//Q;) < -20 A, indicated as thick 
dots (in red in the online version), corresponding to a total 
of 694 stars. Note that, as customary, negative values of the 
equivalent width indicate emission lines. Since at tempera- 
tures Teff > lOOOOK or colours V-I <Q the sample could be 
contaminated by Be stars that are evolving off the MS, there 
we set a more stringent condition on the equivalent width, 
namely Weq(//a) < -50 A. This limit is suggested by a sur- 
vey of the Ha equivalent width of about 100 Be stars in the 
Galaxy (Cote & Waters 1987) in which only one star is found 
with Weq(//a) < -50 A and the largest majority have values in 
the range from -4 A to -30 A. In fact, as we discuss in Sec- 
tionlH all the objects with T^ff > lOOOOK and an Ha excess 
above 4 a are most likely PMS objects with masses around 
3Mq. 

Although we do not have spectroscopic observations of our 
stars to directly confirm our photometrically inferred equiva- 
lent widths, similar studies in the literature have demonstrated 
the success of this approach. White & Ghez (2001) carried 
out a detailed study of 44 T Tauri stars in binary systems in 
Taurus-Auriga, using both photometry and spectroscopy with 
the HST. Their analysis shows that photometric and spectro- 
scopic measurements of the equivalent width of the Ha emis- 
sion line in these objects are in good agreement with one an- 
other More recently, Barentsen et al. (2011) carried out a 
Ha survey of T Tauri stars in IC 1396, a region with consid- 
erable nebular emission. Their photometric determination of 
Weq(//a) for 109 objects in common with the spectroscopic 
survey of Sicilia-Aguilar et al. (2005) also shows an excel- 
lent agreement, particularly for stars with Weq(//a) < -10 A, 
which are the objects of interest in our study. 
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Figure 5. Colour-magnitude diagram of the field of NGC 346 showing as 
open circles the 791 PMS star candidates with excess emission in Ha at the 
4(T level or higher. A total of 694 of them (shown as filled circles, in red in 
the online version) also have Wcq(//a) < -20 A or < —50 A for stars hotter 
than 10000 K. Squares correspond to objects with Wcq(Ha) > -50 A that are 
potential Be stars. The dashed line coiTesponds to the 50 % completeness 
limit of the photometry (see Sabbi et al. 2008). 

4. PHYSICAL PARAMETERS OF THE PRE-MAIN SEQUENCE STARS 

When the bona-fide PMS objects identified in the previous 
section are placed on the CMD (see Figure|5]l, they clearly 
reveal two distinct groups made up of a similar number of 
objects. One is located above (i.e. brighter and redder than) 
the MS and thus is indicative of a population of young PMS 
objects. This population was originally noticed by Nota et al. 
(2006) and later studied by Gouliermis et al. (2007) and Sabbi 
et al. (2007), although these earlier works could only iden- 
tify candidate young (< 5 Myr) PMS stars since no informa- 
tion was available to verify and characterise the PMS nature 
of each individual object. This is the reason why the works 
above could not identify the second group of bona-fide PMS 
objects, which overlap in the CMD with older MS stars in the 
field of the SMC. Hennekemper et al. (2008) detected some of 
these objects using a preliminary version of the colour trans- 
formations between V, I and R that we developed for Paper! 
to search for stars with R-Ha excess in this field, but they did 
not realise what was the nature of these objects. As we show 
here below, these objects are older PMS stars with an age in 
excess of 8 Myr, of the same type as those that we discov- 
ered in the field around SN 1987A in the LMC (see Paper!). 

The physical parameters of the bona-fide PMS stars that we 
have identified are obtained as explained in the following. 

4. 1 . Masses and ages 

The radius comes from the luminosity and effective tem- 
perature of the stars, which in turn follow from the observed 
colour (y - /) and magnitude (V), properly corrected for in- 
terstellar extinction as explained in Section lZTI The stellar 
mass and age were derived by comparing the location of 
each star in the Hertzsprung-Russell (H-R) diagram of Fig- 
ure|6] with the PMS evolutionary tracks and corresponding 
isochrones, via interpolation over a finer grid than the one 
shown in the figure. We followed the interpolation proce- 
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Figure 6. Hertzsprung-Russel diagram of the PMS candidate stars. All 
objects shown here have excess emission in Ha at the 4cr level or higher 
Those indicated with filled circles (in red in the online version) also have 
Wcq(//Q!) < -20 A or < -50 A for stars hotter than 10 000 K. Squares cor- 
respond to objects with Wcq(//o) > —50 A and, as such, are potential Be 
stars. Thick solid lines show the evolutionary tracks from the Pisa group 
(DeglTnnocenti et al. 2008; TogneUi et al. 2011) for metallicity Z = 0.002 
and masses from 0.5 to 4Mq , as indicated. The corresponding isochrones 
are shown as thin dashed lines, for ages of 0.125, 0.25, 0.5, 1, 2, 4, 8, 16 
and 32 Myr from right to left. Note that the constant logarithmic age step has 
been selected in such a way that the typical distance between isochrones is 
larger than the photometric uncertainties. 

dure developed by Romaniello (1998), which does not make 
assumptions on the properties of the population, such as the 
functional form of the IMF. On the basis of the measurement 
errors, and taking into account its position in the H-R dia- 
gram, this procedure provides the probability distribution for 
each individual star to have a given value of the mass and age 
(the method is conceptually identical to the one presented re- 
cendy by Da Rio et al. 2010). 

As for the theoretical models, we adopted those of the 
Pisa group (Degl'Innocenti et al. 2008; Tognelli et al. 
2011), for metallicity Z = 0.002 or about Z0/8. These 
new PMS tracks were specifically computed for the mass 
range from 0.45 M0 to 5 M0 with an updated version of the 
FRANEC evolutionary code (see Chieffi & Straniero 1989 
and DeglTnnocenti et al. 2008 for details and Cignoni et al. 
2009 for an application). 

As mentioned in the Introduction, we opted for a lower 
value of the metallicity, namely Z = 0.002 in place of the 
more customarily adopted Z = 0.004, because the correspond- 
ing tracks provide a better match to the observations, partic- 
ularly where they reach the MS. As previously mentioned, 
Z = 0.002 is still compatible with the metallicity range found 
in the literature for the SMC and is in any case better suited 
to NGC 346 than the value adopted in previous studies of its 
PMS population that made use of the evolutionary tracks of 
Siess etal. (2000) for Z = 0.01. 

Figure|6]reveals that the masses of the bona-fine PMS stars 
span over a decade, from < 0.45 M0 for the coolest objects 
to ^ 4M0 for the hottest ones. As mentioned above, some of 
the objects with T^s > lOOOOK and Weq(//a) > -50 A could 
be Be stars and we have marked all objects of this type with 
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Figure 7. Histograms of the physical parameters of the bona-fide PMS stars 
in the field of NGC 346, namely mass (a), age (b). Ha luminosity (c) and 
mass accretion rate (d). The solid line refers to the 680 objects that also 
have an equivalent width Wc^iHa) < -20 A (or < -50 A for stars with T^ff > 
10000 K) and a well defined mass. For comparison, the dotted line in panel 
(c) refers to the 791 stars with an Ha excess at the 4cr level but without any 
further restrictions on mass, age or equivalent width. 

squares to distinguish them from the rest. Figure|6]also pro- 
vides a more quantitative characterisation of the age differ- 
ence between the two populations of PMS stars that we iden- 
tified earlier in the CMD (Figure|5). The isochrones in the 
figure, corresponding from right to left to ages of 0. 125, 0.25, 
0.5, 1, 2, 4, 8, 16 and 32 Myr, show that the young population 
has a median age of ^ 1 Myr, whereas the older bona-fide 
PMS stars span a wider age range from ^ 10 Myr to 30 Myr, 
with a broad peak around 20 Myr 

We derived a reliable measure of the mass and age for all 
stars indicated as thick dots in Figure|6] except for the coolest 
objects that would have required a more uncertain extrapola- 
tion beyond the limits of the theoretical models. Histograms 
with the mass and age distribution for the 680 stars measured 
in this way are shown in Figures]?^ and |7j), respectively. As 
Figure|6] already implied, the histogram in Figure|7j) shows 
that there is a remarkable paucity of stars with ages around 4- 
8 Myr In particular, the solid line in Figure|7j) provides a his- 
togram of the age distribution with constant logarithmic step 
(a factor of 2) and shows two peaks, at ^ 1 and 20 - 30 Myr 
The two peaks remain visible also when the number of stars 
in each bin is rescaled by the size of the bin, as indicated by 
the dot-dashed line (note that at the extremes of the distribu- 
tion it becomes more difficult to assign an age to the stars, 
so the first and last bin are drawn with a dotted line to indi- 
cate a larger uncertainty). In a companion paper (De Marchi, 
Paresce & Sabbi 201 1), we discuss the age and spatial distri- 
butions of these stars in more detail and show that between 30 
and 10 Myr ago a large number of PMS stars formed, similar 
to that produced by the current burst, which however appears 
to be much more intense. 

4.2. Mass accretion rates 

Besides the Ha luminosity L(Ha) already discussed in the 
previous section, the other physical parameters of interest in 
this work are the accretion luminosity L^cc and the mass ac- 
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cretion rate Mace, that we derived following the procedure 
described in Paper! From an analysis of the literature val- 
ues of Lace and L(Ha) measurements as recently summarised 
by Dahm (2008) for PMS stars in the Taurus-Auriga asso- 
ciation, that paper concluded that the ratio L.^^^/L{Ha) can 
be assumed to be constant. The proportionality constant, ob- 
tained from an elementary fit to the data in the compilation of 
Dahm (2008), is: 



logLaee = ( 1 .72 ± 0.47) + logL(i/a) 



(4) 



where the large uncertainty on the proportionality factor 
arises because the two quantities were determined from non- 
simultaneous observations. Considering that the intensity of 
the Ha line from PMS sequence stars is known to vary by 
about 20 % on a time scale of a few hours and by as much as a 
factor of 2 - 3 in a few days (e.g. Fernandez et al. 1995; Smith 
et al. 1999; Alencar et al. 2001), the large spread can be un- 
derstood. However, as we will show in Section|5] the small 
dispersion of the mass accretion rate with stellar age suggests 
an uncertainty smaller than 0.25. 

The mass accretion rate is related to Laee via the free-fall 
equation, that links the luminosity released in the impact of 
the accretion flow with the rate of mass accretion Mace, ac- 
cording to the relationship: 



Lac 



GM^Mac 



(5) 



where G is the gravitational constant, M, and the mass and 
photospheric radius of the star as derived above and Rm the 
inner radius of the accretion disc. The value of is rather 
uncertain and depends on how exactly the accretion disc is 
coupled with the magnetic field of the star. Following Gull- 
bring et al. (1998), we adopt Ri^ = 5R^ for afl PMS objects. 
By substituting Equation|4] in Equation|5] the mass accretion 
rate Maee, in units of Mq yr~', can be directly expressed as a 
function of L(Ha): 



log 



Ma, 



Mgyr-i 



-7.39 + log ^ + log - log ^ 

^Lq ^Rq ^Mq 



(6) 



L(Ha) R^ M* 
= (-5 .67 ± 0.47) + log + log ^ - log tt^ 

Lq Rq Mq 



The histograms showing the distribution of L(Ha) and Maee 
as derived with Equation|7] for the 680 bona-fide PMS stars 
with well defined masses are shown in Figures|7}; and|7]l, re- 
spectively. The median values of L{Ha) and Maee are, respec- 
tively, 4.6 x 10^' erg/s (or 0.012 Lq) and 3.9 x 10-** Mq yr"'. 
The latter value is about 50 % higher than the median mass ac- 
cretion rate of a population of 133 PMS stars that we detected 
in the field of SN 1987 A (see Paper!) and this is fully consis- 
tent with the much younger average age of the PMS stars in 
NGC346. 

Paper! provides an extensive discussion of the statistical 
and systematic uncertainties involved in determining Mace 
with this method and we refer the reader to that work for more 
details. The typical combined statistical uncertainty on Mace 
for the NGC 346 field is 13 %. !t is largely dominated by the 
uncertainty on L(Ha) of 10 % which in turn is mostly due to 
random errors in the photometry (we note that this uncertainty 
is lower than for the study of the field around SN 1987 A in 
Paper! because of the higher quahty of the ACS photometry). 



As for systematic effects, those on M, and reddening add 
up to 15 %, but it is the rather uncertain relationship linking 
L{Ha) to the accretion luminosity Laee that plays a dominant 
role. We stress here that our method adopts the calibration 
provided by Dahm (2008), which is based on a spectroscopic 
study of Galactic PMS stars. !f the uncertainty on the rela- 
tionship between L{Ha) and Laee is as large as Dahm's (2008) 
work suggests (but see Section|5]l, it could possibly cause a 
systematic effect of up to a factor of 3 on the value of Mace 
that we obtain. 

This is the unfortunate consequence of our incomplete un- 
derstanding of the physics of the accretion process. Improv- 
ing this situation would require direct measurements of the 
bolometric accretion luminosity Laee and simultaneous mea- 
surements of the hot continuum excess emission and of the 
emission line fluxes in Ha or other lines (e.g. Pa/? and Br7). 
!t is also true, however, that these systematic effects plague 
in the same way all determinations of Maee that make use of 
emission lines diagnostics. This includes detailed studies of 
the profile and intensity of the Ha line (e.g. Muzerolle et al. 
1998), since they rely on the same measurements for their cal- 
ibration. Therefore, it is possible and meaningful to compare 
the values of Mace that we obtain here to one another and to 
those derived with other methods, as we do in the following 
section. 

5. EVOLUTION OF THE MASS ACCRETION RATE 

!n this section we investigate how Maee varies with stellar 
parameters such as stellar mass and age. !t is important to 
underline that this analysis is based only on stars that, at the 
time of the observations, were undergoing active mass accre- 
tion as revealed by Ha excess emission above our detection 
threshold. Actually, PMS objects are known to display large 
variations in their Ha emission over hours or days (e.g. Fer- 
nandez et al. 1995; Smith et al. 1999; Alencar et al. 2001), 
with only about one third of them at any given time being 
active Ha emitters with WeqiH) < -10 A in the field around 
SN 1987A (Panagia et al. 2000). 

!n fact, a similar requirement of a high accretion activ- 
ity is common to all studies in which PMS stars are identi- 
fied based on single-epoch observations, including those con- 
ducted spectroscopically on T Tauri objects in the Milky Way 
(e.g. Sicilia-Aguilar et al. 2006). This implies that the ac- 
cretion rates measured for the detected objects will, in fact, 
represent the highest values possible for that class of objects, 
while many more stars may well be accreting at more modest 
rates. As a result, the true average accretion rate appropriate 
for a given class of objects will be just a fraction of the values 
measured for the members of that class that can be detected. 
On the other hand, since all studies made with different tech- 
niques suffer from the same kind of selection on the accretion 
efficiency, a direct comparison between their results is still 
possible and meaningful. 

5.1. A qualitative analysis 

Let's start with a quaUtative analysis of the H-R diagram 
of the bona-fide PMS stars shown in Figure|8^, where the size 
of the symbols is set to be proportional to Maee according to 
the legend. Besides the 680 bona-fide PMS stars in NGC 346 
(darker dots, red in the online version) we have included in 
the figure for comparison also the 133 bona-fide PMS stars 
in the field of SN 1987 A discussed in Paper! (lighter dots, 
gray in the online version). The theoretical PMS isochrones 
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Figure 8. H-R diagrams of all 680 bona-fide PMS stars. In panels (a), (h) 
and (c) the size of the symbols is set to be proportional respectively to Mace , 
L(Ha) and Weq(//cj), according to the legends. The darker symbols (in red 
in the online version) correspond to NGC 346 (same as Figure|6), while the 
lighter ones (gray in the online version) represent the bona-fide PMS stars 
around SN 1 987 A. The isochrones corresponds from right to left to ages from 
0.125 Myr to 32Myr with a constant logarithmic step and are for metallicity 
Z = 0.002. Therefore, isochrones are only applicable to the NGC 346 data. 



are the same as shown in Figure|6]and correspond to a metal- 
Ucity Z = 0.002 and to ages from 0.125 Myr to 32 Myr in- 
creasing from the right with a constant logarithmic step of a 
factor of 2. While not directly applicable to the PMS stars 
in the SN 1987 A field, due to the different metallicity, these 
isochrones serve nonetheless the purpose to illustrate the ef- 
fects of increasing PMS age. 

It is noteworthy that the mass accretion rate varies by over 
three orders of magnitude across the H-R diagram in a very 
regular way: younger stars, which are farther away from the 
MS, have typically higher mass accretion rates, as expected 
from the theory of stellar evolution. Also the stellar mass ap- 
pears to affect the value of Mace, in the sense that less massive 
stars tend to have smaller accretion rates. It is also quite re- 
markable that even PMS stars very close to the location of the 
zero-age MS (ZAMS), which for the masses of interest here 
can be considered equivalent to the 32 Myr PMS isochrone, 
are still undergoing substantial mass accretion. 

To verify that this is not a spurious effect due to errors in 
our derivation of Mace, we show in Figure[8]3 the same H-R 
diagram but with the size of the symbols proportional to the 
measured Ha luminosity. There one sees that L(Ha) can be 
rather large for objects at or near the ZAMS, suggesting that 
the accretion process remains relatively active when the star 
approaches the MS. More precisely, this demonstrates that, 
even if Mace decreases as a PMS star approaches the MS, its 
shrinking radius allows the star to extract more potential en- 
ergy per unit rate of mass accretion (i.e. for a given rate in 
gem"'), thereby making the object easier to detect in Ha. Fi- 
nally, Figure|8}; shows that the EW{Ha) values of PMS ob- 
jects near or at the MS are typically much larger than those 
caused by chromospheric activity for stars in this mass range 
(see e.g. Ulmschneider 1979; Pasquini et al. 2000), confirm- 
ing that they are bona-fide young objects. 

5.2. A quantitative analysis 
For a more quantitative analysis of the evolution of the mass 
accretion rate, we turn to Figure|9]where the variation of Mace 
is shown as a function of stellar age. All 680 bona-fide PMS 
stars in NGC 346 are indicated as diamonds, with a symbol 
size proportional to their mass. Besides being very homo- 
geneous, this sample is about five times larger than all Mace 
measurements today available for Galactic PMS objects taken 
together (e.g. Sicilia-Aguilar et al. 2006). Furthermore, it 
spans over a decade in mass and two orders of magnitudes in 
age while including objects from the same star forming region 
with very accurate and uniform measurements. All this makes 
the sample particularly well suited for a quantitative analysis 
of the dependence of the mass accretion rate on stellar param- 
eters. 

The dashed line in Figure|9]represents the best fit to the evo- 
lution of the mass accretion rate in NGC 346, with a slope 
a = -0.54 and intercept Q = -l (at 1 Myr). The fit appears 
considerably less steep than the prediction of Hartmann et al. 
(1998; see also Muzerolle et al. 2000), based on models of 
viscous disc evolution and represented here by the solid line, 
which however seems to adequately reproduce the trend of de- 
creasing Mace with stellar age for low-mass Galactic T-Tauri 
stars as compiled by Sicilia-Aguilar et al. (2006) and also 
shown in the figure. The different slopes of the observed and 
model relationships translate into predicted Mace values that 
are about two order of magnitudes smaller than what we ac- 
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Figure 9. Mass accretion rate as a function of stellar age. The 680 bona-fide 
PMS stars that we identified in NGC 346 are indicated as diamonds, with a 
symbol size proportional to the stellar mass. The dashed line is the best fit to 
their distribution, with slope a = -0.54 and intercept Q = -l (at 1 Myr). The 
other symbols (see legend) correspond to the Galactic T-Tauri stars studied by 
Sicilia-Aguilar et al. (2006), whose distribution shows a lai'ge scatter but is 
nonetheless consistent with the models of viscous disc evolution of Hartmann 
et al. (1998) represented by the solid line. 

tually observe for an age of 10 Myr. 

This discrepancy had already been noticed in the field of 
SN 1987 A studied in Paper I, although it is pointed out there 
that the Mace value measured at 10 Myr is in excellent agree- 
ment with that of the most massive PMS stars (with spec- 
tral type G) in the Trumpler37 sample of Sicilia-Aguilar et 
al. (2006) shown as large circles in the figure. Even though 
the age of the G-type stars in Trumpler37 is very uncertain 
(Sicilia-Aguilar et al. 2006) and could be much younger (J. 
MuzeroUe, priv. comm.; see also Hartmann 2003), we can- 
not exclude that the slower decline of Mace with time that we 
observe in NGC 346 and around SN 1987 A could be due to 
the typically higher mass of our objects compared with that 
of the Galactic T Tauri stars that were the subject of recent 
studies. A dependence of Maee on the stellar mass M is in- 
deed expected and could be as steep as Mace M^ according 
to MuzeroUe et al. (2003; 2005), Natta et al. (2004; 2006) 
and Calvet et al. (2004), although the uncertainties are large 
(see also Clarke & Pringle 2006 for a different interpretation). 

Owing to the large size and wide mass range of the objects 
in our sample, we can investigate the mass dependence of 
Mace in a very robust way. As a first step, we divide out sam- 
ple in four mass groups, namely O.4-O.8M0 , 0.8- 1.1 M0 , 
l.l-1.5M0and 1.5-3M0, and plot for each one sepa- 
rately the temporal evolution of Mace in the four panels of 
Figure[TO] The slope a and intercept Q (at 1 Myr) of the best 
linear fit to the data, according to the relationship log Mace = 
a X log(age) + Q, are given in each panel together with their 
uncertainty. It appears that the slopes of the four mass groups 
are in very good agreement with one another, within the un- 
certainty of the fits, but remain significantly different from the 
trend predicted by Hartmann et al. (1998) and MuzeroUe et 
al. (2000) for viscous disc evolution (solid lines). 

In principle, the shallower slopes of our fits could be af- 
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Figure 10. Same as Figure|9]but for stars in four mass groups, as indicated 
in each panel. We also provide the slope a and intercept Q (at 1 Myr) of 
the best fits (thick long-dashed lines). The thin short-dashed lines correspond 
to our L(Ha) detection limits, which appear to not significantly affect the 
determination of the slopes, except possibly for the least massive stars at ages 
> 10 Myr 

fected by our very strict criteria for selecting bona-fide PMS 
stars. Requiring that the Ha excess be at the 4cr level or 
higher and that Weq < -20 A will set a lower limit to the Ha 
luminosity that we can accept as significant and, therefore, 
to the smallest Maee value that we can detect at a given mass 
or age. If this lower limit is too high, we could be ignoring 
too many objects with low mass accretion rate, particularly at 
older ages, and therefore skew the measurement of the slope 
towards low values. On the other hand, it appears that our 
detection limits (shown by the thin short-dashed lines for the 
average mass of each group) do not significantly affect the de- 
termination of the slope, except possibly for the least massive 
stars at ages > 10 Myr For this group, however, the slope is 
determined by the much more numerous young stars. 

More importantly, the upper envelopes of the distribution 
of all four mass groups appear to be fully consistent with the 
slope of the best fit. This is shown graphically by the thin 
long-dashed lines that represent the best-fitting line shifted 
vertically by 0.5 dex. While our conservative limits on the Ha 
excess emission may prevent us from detecting objects with 
low L(Ha), no selection effects are to be expected for stars 
with high Ha luminosity. Since the distribution of these ob- 
jects is consistent with the slope of the best fit, the latter must 
be correct. Therefore, within the uncertainties quoted in Fig- 
urefTOl we can conclude that all four mass groups are compati- 
ble with the same slope a ~ -0.6 (formally a = -0.55 ±0.02). 

An interesting feature of Figure[TO] is that the spread of 
the logMacc data values around the best fitting lines is small, 
typically ^ 0.25 dex (rms). This proves quite convincingly 
that the true uncertainty in the relationship between log Lace 
and log L(//a) given by Equation]?] cannot be as large as the 
0.47 dex value suggested by our elementary fit to the data in 
the compilation of Dahm (2008), but it has to be intrinsically 
smaller than the ^ 0.25 dex factor that we observe in our own 
NGC 346 data. We conservatively assume 0.25 dex (or a fac- 
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Figure 11. Run of Macc/"i as a function of age for the PMS objects in 
NGC 346 (circles) and for Galactic stars in Taurus (crosses) and in Trum- 
pler 37 (squares). The hatched band show the best fit (±1 ct) to the NGC 346 
data, the dark-shaded band is the fit to the PMS stars in the 30 Dor region as 
derived in Paper III and the light-shaded band is the fit to the Galactic stars 
(including also the Orion Nebula Cluster; Da Rio et al. in preparation). 

tor of < 1 .8 in Hnear units). 

Having estabhshed that the data in all four panels of Fig- 
urefTOl are consistent with the same slope a ~ -0.6, we can 
use the intercept 2 at 1 Myr to study how Mace varies with the 
stellar mass m. It appears that the value of Q in the four panels 
of Figure[TO]decreases with increasing mass approximately as 
m"' . The formal uncertainty on Q may seem too large to reach 
this conclusion, but the scatter is due to the grouping of stars 
of different mass in the same bin. It is, therefore, more ap- 
propriate to perform a multivariate fit to the three variables 
simultaneously (Mace, age and mass). We opt for a linear mul- 
tivariate least-square fit of the type: 

logMace = fl X log f + X log m + c (7) 

where f is in years, m the mass in solar units and c is a 
constant. The formal best fit gives a = -0.61 ±0.02, b = 
0.90 ± 0.10 and c = -3.27 ± 0.24. However, since the set of 
parameters a = -0.6, b = 1 and c = -3.3 would give almost 
equally small residuals, we will take them hereafter as a good 
as well as convenient representation of the dependence of Mace 
on stellar mass and age in NGC 346. 

5.3. Comparison with other star fanning regions 
From this analysis it appears that the NGC 346 data are not 
consistent with a steep dependence of Mace on the stellar mass 
of the type reported by Muzerolle et al. (2003) and Calvet 
et al. (2004) in nearby star forming regions. However, the 
gentler decline of Mace with mass that we find here is fully 
compatible with that obtained in a study of about 900 PMS 
stars recently identified in the regions around 30 Dor (Spezzi 
et al., in preparation. Paper 111) as well as with the preliminary 
results of a ground-based study of about 500 PMS objects in 
the Orion Nebula Cluster (Da Rio et al., in preparation). In all 
these cases, values of a = -0.6 and b = I in Equation|7]give a 
good fit to the data. 



This suggests that we could use the same functional form 
and parameter values to fit the observed evolution of Mace for 
Galactic PMS stars and compare the result to the data in the 
Magellanic Clouds. We do so in Figure[TTl where we show the 
run of the specific mass accretion rate, defined as m = M^^c/ni, 
as a function of time for the PMS objects in NGC 346 (cir- 
cles) and for Galactic stars in Taurus (crosses; from Calvet 
et al. 2000) and Trumpler 37 (squares; from Sicilia-Aguilar 
et al. 2006). The hatched band represents the best fit to the 
NGC 346 data, the dark-shaded band is the fit to the PMS stars 
in the 30 Dor region as derived in Paper 111, whereas the light- 
shaded band refers to Galactic stars (including also the Orion 
Nebula Cluster; Panagia et al. in preparation). The hatched 
and solid bands correspond to the ±1 cr scatter (i.e. the range 
from 17 to 83 percentile) and indicate that the SMC and LMC 
data are compatible with one another, but that the difference 
with Galactic PMS stars is significant. 

The logarithmic value of the specific mass accretion rate at 
1 Myr, q = logm(l Myr) calculated using the best fitting line 
as given in Equation|7] can serve as a measure of this dif- 
ference, q takes on the values of -6.9 ± 0.2, -6.7 ± 0.4 and 
-7.7 ±0.6 respectively for NGC 346, the 30 Dor region and 
the Galaxy. As akeady pointed out in Paper 1, the higher Mace 
values of PMS stars in the LMC compared to Galactic ob- 
jects of the same age could result from the lower metallicity 
of the clouds. The radiation pressure of the forming star is 
expected to be less strong on lower-metallicity disc material 
and this will delay the dissipation of the disc, thereby keeping 
the accretion process active for a longer time. This scenario is 
compatible with the PMS stars in the Galaxy having a smaller 
value of the specific mass accretion rate q than those in the 
Magellanic Clouds. 

A possibility would be that the shallower photometric depth 
of the LMC data may cause the L{Ha) detection threshold for 
PMS stars to be higher (see Paper 111 for more details) and 
thus could skew the median value of Mace- We will address 
this possibility in detail through the analysis of recently ob- 
tained HST observations of the same region (De Marchi et 
al., in preparation). Alternatively, it is also possible that be- 
low a certain metallicity threshold the radiation pressure on 
the disc material no longer changes. Investigating this pos- 
sibility requires the analysis of larger and homogeneous data 
set of Mace measurements in low-metallicity environments. 

5.4. Mass growth during the PMS phase 

Finally, we would like to highlight an interesting corol- 
lary result stemming from the large mass accretion rate val- 
ues that we find for the objects in NGC 346. Using the 
rates as expressed by Equation|7] we can estimate the total 
mass accreted by a star over the time span from birth to the 
ZAMS. Already in 1993 Palla & Stahler (1993) showed that 
the timescale to reach the ZAMS is a strong function of the 
stellar mass. Fitting the more modem evolutionary calcu- 
lations (DeglTnnocenti et al 2008; Tognelli et al 2011) for 
metallicities I/2-I/8Z0 and stellar masses in the range 0.6 
up to 4 Mq , to within 8 % accuracy we can express that time 
in years as: 

fzAMS =i 20 X 10^ (Z/0.002)"-5 m-^ l (8) 

Equation|2] provides an analytical expression for the aver- 
age mass accretion rate of the detected stars. On the other 
hand, there must be an appreciable fraction of PMS objects 
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that, at any given time, are not detectable as strong Ha emit- 
ters. Therefore, in order to properly estimate the amount of 
accreted mass over the PMS phase we have to make allowance 
for the appropriate duty cycle. While doing this would require 
a detailed knowledge of the time evolution of the Ha emission 
for a large sample of PMS stars in this region, for the sake 
of simplicity we will assume that the evolution consists in a 
series of two-level transitions that recur in time. In this sim- 
plified model the high state value corresponds to the average 
value measured for the detected PMS stars and the low state 
of the Ha undetected stars is zero. Thus, the true mass accre- 
tion rate of any star of mass m at the time t will simply be the 
average value measured for the detected PMS stars multiplied 
by an efficiency factor that in principle is a function of the 
stellar mass, the time and the metallicity of the region. 

An approximate value of this function is given by the ra- 
tio of the number of PMS stars with detected Ha emission 
and the total number of candidate PMS stars as defined by 
their location in the H-R diagram. In a companion paper, De 
Marchi, Panagia & Sabbi (2011) show that for NGC 346 such 
a ratio is essentially constant over time up to PMS ages of 
8Myr at a level of < >~ 0.28. Note that similar values 
have been found in the field of SN 1987 A for PMS stars with 
median age of 14Myr (< </> >~ 0.32; Panagia et al 2000). 
Therefore, adopting a constant value of < >= 0.28 for PMS 
stars in NGC 346 allows us to integrate the accretion rate over 
time from up to fzAMS in order to estimate the total amount 
of mass (in M0 units) accreted by these objects during their 
PMS phase. From Equations|7]and 8 we can express it in units 
of M0 as follows; 

Macc= I <(/» Macct/? (9) 

rfZAMS 

= / <</» t" 10" dt 
Jq 

= 4x 10""^ «t» m / r^^-^dt 
Jo 

= 4xlO-4<0>m^ 



= 0.23 m« »^ 

showing that the amount of accreted mass is virtually insen- 
sitive to the stellar mass. In the range 0.4 M© < m < 4 Mq as 
studied here, EquationfTOl implies M^cc = 0.23 ±O.O3M0 . 

Taken at face value, such a high Mace figure would imply 
that a considerable amount of mass is accreted during the 
PMS phase, at least in low-metallicity environments such as 
the Magellanic Clouds. If this is the case, the PMS evolu- 
tion of moderate-mass stars, say < 2 Mq , should be reconsid- 
ered and recalculated taking into account the high Mace val- 
ues. One of the effects of Equation 9 is that the appropriate 
fzAMS for a star that reaches the MS with a given mass will be 
longer than the value estimated by evolutionary models that 
assume Maee = after the first few 10"^ yr (e.g. Baraffe et al. 
2009). Since the star had an appreciably lower mass for a very 
long time (several Myr) at the early epochs, it was necessarily 
evolving more slowly. 

Extrapolating these findings to lower masses than those that 
we have reached here (i.e. < 0.4 M© ) would carry the impli- 



cation that no star could form with a mass below a suitable 
minimum value. For NGC 346 this would correspond to a 
lower cut-off of the initial mass function at ~ 0.2 M0 . In 
the Milky Way, where Maee is about five times lower but the 
metallicity is about 8 times higher so that the PMS evolution- 
ary times are a factor of about 2.8 longer, the cut-off would 
occur at ^ 0.06 M© . 

On the other hand, it is also possible that for low-mass stars 
circumstellar discs are eroded in a finite time that is apprecia- 
bly shorter than their fzAMS- For example, if for stars below 
1 M0 we were to limit the integrals in Equation 9 between 
and the smaller of fzAMS and 20 Myr, the accreted mass would 
decrease proportionally to the mass, i.e. 

Maee ^ 0.23 OT for m<lM0. (10) 

In practice this would still require some refinements to 
models of stellar evolution in the PMS phase, but the increase 
in mass would never be higher than 20 % of the final ZAMS 
mass. Note that the choice of 20 Myr as the upper limit for 
the integration stems from the fact that near the MS we still 
find 1 Mq stars that have strong Ha excess, thus indicating 
the presence of considerable accretion. Deeper observations 
of NGC 346 reaching well below 1 Mq on the MS can eas- 
ily show whether the discs of lower-mass objects are indeed 
eroded before they reach the ZAMS. It would be sufficient 
to verify whether strong Ha excess emission is still shown by 
stars smaller than 1 Mq near the MS, or, alternatively, whether 
for objects of these lower masses strong Ha emission is only 
seen at younger ages, when they are still well above the MS. 
If at these masses the discs are indeed eroded before reaching 
the ZAMS, Ha emission should systematically taper off for 
objects older than a cut-off age. In any case, these new results 
can have important implications for some of the assumptions 
and boundary conditions used in theories of star formation 
and early evolution. 

6. SUMMARY AND CONCLUSIONS 

We have studied the properties of the stellar populations 
in the field of the NGC 346 cluster in the SMC as obsei-ved 
with the ACS camera on board the Hubble Space Telescope 
(Sabbi et al. 2007). To this aim, we have employed a novel, 
self-consistent method that allows us to reliably identify PMS 
stars undergoing active mass accretion, regardless of their age. 
The method, fully described in Paper 1, does not require spec- 
troscopy and combines broad-band V and / photometry with 
narrow-band Ha imaging to identify all stars with excess Ha 
emission and allows us to derive the accretion luminosity Laee 
and mass accretion rate Maee for all of them. The main results 
of this work can be summarised as follows. 

1. From the photometry of Sabbi et al. (2007), we have 
selected all objects in a region of 200" x 200" around 
the centre of NGC 346 that have a combined mean error 
6i in the three bands (V, I and Ha) of less than 0.1 mag, 
as defined in Equation[T] A total of 18 764 stars satisfy 
this condition. 

2. The CMD reveals the presence of both young and old 
populations, but the conspicuous broadening of the up- 
per MS proves that there is also variable extinction. 
We have determined the reddening towards each of 
these stars from their colour and magnitude displace- 
ments with respect to a 4 Myr isochrone for metallicity 
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Z = 0.002 and a distance modulus of 18.92. We found 
that extinction is in the range 0.1 < E(V-I) < 0.25 or 
0.27 <Av < 0.67 (respectively 17 % and 83 % Hmits), 
with a median value of E(V -I) = 0.16 or Ay = 0.43. 
We have then used the reddening values towards each 
of the selected stars to derive a reddening correction for 
all objects in their vicinity. 

3. Following the method developed in Paper I, we have 
identified 791 PMS candidates as being objects with 
Ha excess above the 4 a level with respect to the ref- 
erence provided by normal stars. Their average Ha lu- 
minosity is 2.7 X lO-'' ergs"' or ^ 10"^ Lq. We have 
also determined the equivalent width of the Ha emis- 
sion line of these objects and have classified as bona- 
fide PMS stars all those with Weq < -20 A (or < -50 A 
for stars with Teff > lOOOOK). These conditions guar- 
antee that our sample is free from contamination due to 
the chromospheric activity of older objects or to the ro- 
tational winds of Be stars. A total of 694 objects satisfy 
these conditions. 

4. By comparing the locations of these objects in the H- 
R diagram with the PMS evolutionary models of the 
Pisa group for metallicity Z = 0.002, we have been 
able to derive an accurate value of the mass and age 
of 680 bona-fide PMS stars (note that previous deter- 
minations of these parameters for candidate PMS stars 
in NGC 346 that made use of evolutionary models for 
Z = 0.01 are necessarily less accurate). The masses of 
these objects range from 0.4 Mq to 4 Mq , with an aver- 
age value of ^ 1 M0 . Their ages show a clear bimodal 
distribution with two peaks at ^ 1 Myr and ^ 20 Myr 
and very few objects around ^ 8 Myr, revealing the 
presence of two distinct but equally populous genera- 
tions of stars. We address in a companion paper (De 
Marchi, Panagia & Sabbi 2011) the properties of and 
the relationships between the two generations. 

5. From the Ha luminosity and the other physical param- 
eters that we have measured, we have derived the mass 
accretion rate M^cc of all bona-fide PMS stars. The me- 
dian value of Mace is 3.9 x 10"^ Mq yr"'. This value is 
about 50 % higher than that measured in Paper 1 for a 
population of 133 PMS stars in the field of SN 1987A, 
owing to the much younger median age of PMS objects 
in NGC 346. 

6. Thanks to the unprecedented size of our PMS sample 
and of its spread in mass and age, we have been able to 
study the evolution of the mass accretion rate as a func- 
tion of stellar parameters. Regardless of the mass of the 
star, our analysis shows that the mass accretion rate de- 
creases with roughly the square root of the age, or about 
three times slower than predicted by current models of 
viscous disc evolution, and that more massive stars have 
systematically higher mass accretion rate in proportion 
to their mass. A multivariate linear regression fit reveals 
that logMacc — -0.61ogf + log m + c, where t is the age 
of the star, m its mass and c a quantity that is higher at 
lower metallicity. 

7. The large mass accretion rates that we find imply that 
a considerable amount of mass is accreted during the 
PMS phase, at least in low-metallicity environments 



such as the Magellanic Clouds. At face value, the ob- 
served rates integrated over the entire PMS phase would 
give a total accreted mass of ^ 0.2 Mq , with negligi- 
ble dependence of the final mass on the ZAMS, thereby 
possibly suggesting a lower cut-off mass in the stellar 
mass function of NGC 346. If, however, the circum- 
stellar discs are eroded on a time scale shorter than 
the PMS lifetime, the total accreted mass will be of 
course smaller but not negligible. In any case, the PMS 
evolution of moderate-mass stars (< 2M0 ) should be 
reconsidered and recalculated taking into account the 
high Mace values, since for a given MS mass the evolu- 
tionary time needed to reach the ZAMS will be longer 
than what is currently estimated by models that assume 
Mace = after the first few lO"* yr. 
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